Free-induction decay and envelope modulations in a narrowed nuclear spin bath 
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We evaluate free-induction decay for the transverse components of a localized electron spin coupled 
to a bath of nuclear spins via the Fermi contact hyperfine interaction. Our perturbative treatment 
is valid for special (narrowed) bath initial conditions and when the Zeeman energy of the electron b 
exceeds the total hyperfine coupling constant A: b > A. Using one unified and systematic method, 
we recover previous results reported at short and long times using different techniques. We find a 
new and unexpected modulation of the free-induction-decay envelope, which is present even for a 
purely isotropic hyperfine interaction without spin echoes and for a single nuclear species. We give 
sub- leading corrections to the decoherence rate, and show that, in general, the decoherence rate has a 
non-monotonic dependence on electron Zeeman splitting, leading to a pronounced maximum. These 
results illustrate the limitations of methods that make use of leading-order effective Hamiltonians and 
re-exponentiation of short-time expansions for a strongly-interacting system with non-Markovian 
(history-dependent) dynamics. 

PACS numbers: 03.65.Yz, 72.25.Rb, 31.30.Gs 



I. INTRODUCTION 

The hyperfine interaction of a quantum- dot- confined 
electron spin with surrounding nuclear spins and the re- 
sulting decoherence of electron-spin states has been a fo- 
cus of research in the last several years because of poten- 
tial applications in spintronics^' 2 -' 3 - and quantum informa- 
tion processings^ Hyperfine-induced electron-spin de- 
coherence is one of the most significant obstacles to viable 
quantum computation with confined electron spins. It is 
therefore of central importance to understand this deco- 
herence mechanism so that schemes can be developed to 
suppress it. 

One of the most promising strategics to suppress 
spin dephasing is to prepare the nuclear-spin system 
in a less-noisy 'narrowed' statei^&SilS Once such a 



state is prepared, it can be maintained over an as- 
tonishingly long time scale, exceeding hours yii since 
spin diffusion processes are highly suppressed near 
confined electron spins Recently, great progress 
has been made in experimentally realizing such state 
narrowing ) 11 : 13 : 14 ' 15 : 16 ' 17 : 18 as well as single-spin readout 
and coherent control ) 19 ' 20 ' 21 ' 22 ' 23 which we expect to lead 
to improved coherence-time measurements in the very 
near future. 

In the absence of rcfocusing pulses, and at time scales 
that are short compared to the relevant time scale for 
the nuclear dipolar interaction, an electron spin interact- 
ing with a narrowed nuclcar-spin environment dephascs 
mainly via virtual flip-flops with the nuclear spins, pro- 
vided a sufficiently strong magnetic field is applied to 
suppress direct spin flips between electron and nuclear 
spins. Under these conditions, the electron-spin dynam- 
ics pass through various stages with a zoo of different de- 
cay laws, obtained by various methods (see, e.g., Fig. 5 of 



Ref. 1241 ): an exact solution for a fully-polarized nuclear 
system and leading-order generalized master equation 
(GME) have both shown a short-time (partial) power- 
law decay ; 7 ' 25 ' 26 an effective-Hamiltonian and short-time- 
cxpansion approach shows that the initial partial decay 
is followed by a quadratic decay shoulder , 27 ' 28 and a 
Born-Markov approximation applied to the same effec- 
tive Hamiltonian shows that the majority of the decay 
is typically exponential in the high-field (perturbative) 
regimei 29 ' 30 ' 31 Finally, an equation-of-motion approach 
has shown a long-time power-law decay to zeroi 32 ' 33 

In this article, we show that each of these results can 
be obtained in a systematic way from a single unified 
approach, by extending the GME introduced in Ref. 
to higher order. In addition to recovering previous re- 
sults at all time scales, we find important qualitatively 
new features, including a modulation of the decay en- 
velope (even for a fully isotropic hyperfine interaction). 
Moreover, we give sub-leading corrections (in the inverse 
electron Zeeman splitting 1/6) to the decoherence rate 
I/T2 calculated previously^ These corrections suggest 
an interesting non-monotonic dependence of I/T2 on 6. 
Neither the envelope modulations, nor the sub-leading 
corrections to I/T2 can be found from dynamics under 
the effective Hamiltonian alone. The results presented 
here therefore show limits to the validity of some previous 
approaches based on high-order expansions of a leading- 
order effective Hamiltonian. 

The rest of this article is organized as follows: In Sec. 
UT1 we introduce the relevant Hamiltonian, initial condi- 
tions, and exact equation of motion (GME) for the elec- 
tron spin. Sec. HHI contains a review of the systematic ex- 
pansion for the electron-spin self-energy (memory kernel) 
in powers of electron- nuclear flip-flops Vg , and gives the 
result up to fourth order in Vg. In Sec. IIVI we evaluate 
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the full non-Markovian spin dynamics for an electron in 
a two-dimensional quantum dot. We recover previously- 
known results found using various other methods and 
present new results for the envelope modulation and cor- 
rections to the exponential decoherence rate I/T2. In 
Sec. [V] we explore the behavior of the fourth-order so- 
lution in the non-perturbative regime and comment on 
the range of validity of this technique and possible ex- 
tensions to higher order. We conclude in Sec. I VII with a 
summary of the results found here and a comparison of 
these results with those that have been presented in the 
literature. Technical details are given in Appendixes lAl 
and[Bj 

II. HAMILTONIAN AND GENERALIZED 
MASTER EQUATION 

We consider the Hamiltonian for a localized electron 
spin-1/2 (with associated spin operator S), interacting 
with a bath of nuclear spins Ik via the Fermi contact 
hyperfme interaction. We allow generally for a Zeeman 
splitting b = cji-IbB of the central spin S and site- (or 
species-) dependent Zeeman splitting 67^ in the bath for 
a nuclear spin Ij. at site k. The Hamiltonian for this 
system is (setting H—l) 

H = bS z +6^ 7fc / fe z + S-h; h = ^A k I k , (1) 

k k 

where the hyperfme coupling constant at site k is given 
by A k = vqA^ |^(rfc)| 2 if the nucleus at site k is of iso- 
topic species j k with associated total hyperfme coupling 
constant A* k , ip(r k ) is the electron envelope wavefunc- 
tion, evaluated at site r k (the position of the fc th nuclear 
spin), and vq is the atomic volume. 

In Eq. |T]) we have neglected the nuclear dipole- 
dipole interaction, which can give rise to additional in- 
ternal dynamics in the nuclear spin system, and con- 
sequent decay of the electron spin i 27 i 30 i 34 i 35 i 36 Dipolc- 
dipolc-induced nuclear spin dynamics are highly sup- 
pressed in the presence of an inhomogeneous quadrupo- 
lar splitting^ 7 - or Knight-field gradient in a small quan- 
tum dot (the "frozen-core" or diffusion-barrier effect, see 
Ref. [H for a review). The relevant decoherence rate 
due purely to the hyperfme interaction is enhanced for 
a small dot (I/T2 ~ 1/N for a quantum dot containing 
N nuclear spins), whereas the dipolc-dipolc- induced nu- 
clear dynamics arc suppressed for a small dot due to the 
frozen-core effect. Thus, there will always be some dot 
size TV where the nuclear dipolar interactions can be ne- 
glected, even up to times that are long compared to the 
electron-spin decoherence time. 

A. Initial conditions 

We choose product-state initial conditions 

p(0)=ps(0)®M0), (2) 



where Pi(S) — T r S(j)P is the reduced density matrix for 
the nuclear (electron) system. Such an initial state can 
be prepared through fast strong pulses applied to the 
electron spin or by allowing an electron to tunnel rapidly 
into a localized orbital^ 

Typically, nothing will be known about the nuclear- 
spin system at the beginning of an experiment, and the 
density matrix pi (0) will be well-characterized by a com- 
pletely random (infinite temperature) mixture. Random- 
ized initial conditions for the nuclear-spin bath result in 
a rapid Gaussian decay of the transverse electron spin in 
the presence of a strong Zeeman splitting fri L 19 ' 25 i 39 ' 40 i 41 
This rapid decay, due to static fluctuations in the initial 
conditions, can be removed by performing a measurement 
of the z-component of the slowly-varying nuclear field 
/i 2 fi There have been several theoretical proposal o 8 i 9 i 10 
to measure the nuclear-spin system into an eigenstate of 
h z and there are now several experiments where similar 
state preparation has been achieved through dynamical 
pumping. 11 i 14 i 15 i 16 i 17 i 18 After preparing the nuclei in an 
eigenstate of the operator h z , the nuclear system will be 
described most generally by an arbitrary mixture of de- 
generate /i^-eigenstates \rii): 

Pi(0) = ^2 Pa K) Kl +^2pij \ n i) (n>j\: ( 3 ) 

where 

h z \m) = h z n \m) v». (4) 

In this paper, we will assume that there is no 'special' 
phase relationship between the different /i z -eigenstates, 
which allows us to approximate pi(0) by the diagonal 
part of Eq. © 

pz(0) « ^2pu\rii) (n^ , (5) 

i 

where, for any particular i, the state |n.;) is given by a 
product of I|-eigenstates: 

\ni) = (S}\l jk mi), (6) 
k 

with 7| \ljk m k) = TO fc \ ljk m k) an d where — Tj k < m l k < 
/,. • 

We will find it convenient to define the average of an 
arbitrary function of .P-eigcnvalues f(m): 

((f(m)))=Y / Pu(n\nmi)\n i ) ) (7) 

i 

where we assume a uniformly polarized nuclear-spin sys- 
tem throughout this article, making the average on the 
right-hand side independent of k. 

B. Generalized master equation 

The von Neumann equation for the density matrix 
p = —i [H, p] = —iLp can be rewritten in the form of 
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the Nakajima-Zwanzig generalized master equationi 42 i 43 
If we introduce a projection superoperator P that pre- 
serves the initial condition Pp(0) = p(0), the Nakajima- 
Zwanzig GME can be written as 



Pp(t) = -i?l?p{t)-i f dt'T{t-t')?p{t'), 
Jo 



E(i) 



-iPLQe- lLQt QLP. 



(8) 
(9) 



Additionally, the projector P must satisfy P 2 = P and 
is typically chosen to preserve all system variables S a : 
TrS a p(t) = TrS a Pp(t). Here, Q is the complement pro- 
jector: Q = 1 — P. For the special case of the Hamilto- 
nian (Eq. ([1])) and the initial condition ([5]), we choose 
the projection superoperator P = p/(0)Trj and find that 
the exact equation of motion for the transverse electron 
spin (S± =S X ± iS y ) is of the form 7 . 

j t (S + ) t = iu n (S+) t -i J dt'X(t-t')(S+) t ,, (10) 
S(f) = Tr [S+T(t)S-pi(0)], (11) 



where 



b + K 



C. Rotating frame 



We define the coherence factor xt in the rotating frame 
oh = 2e^ l(w " +Aa,)t (S+) t , (12) 



and the associated self-energy 

E(t) = e -*( w »+ Au )*E(t), 



(13) 



with Lamb shift Aw due to virtual excitations of the 
bath: 



(14) 



Aw = -Re / dt£(t). 



This gives an equation of motion for xt- 



xt = —iAujxt — i / dt'T,(t — t')xf. (15) 
Jo 

Eq. (JT5J) is an exact equation of motion for the coherence 
factor xt, and therefore serves as an important starting 
point for sytematic approximations in the rest of this 
article. 



III. SELF-ENERGY EXPANSION 

In the absence of an exact closed-form expression for 
E(i), we must resort to an approximation scheme. For 
a large electron-spin Zeeman splitting 6, and due to the 
large difference between the magnetic moments of elec- 
tron and nuclear spin (7^ ~ 10 -3 ), it is appropriate to 



separate the Hamiltonian (Eq. (p}) into an unperturbed 
piece that preserves S z and a flip-flop term, which in- 
duces energy non-conserving flip-flops between electron 
and nuclear spins: H = Hq + Vg , where 

H = loS z + b^lklk; u = b + h z , (16) 

k 

1 



= ^{h+S- + h-S+ 



(17) 



We can then write E(t) in powers of Vg by performing 
a Dyson-series expansion of Eq. (0) and inserting the 
result into the definition (Eq. (fTTj) );! 



(18) 



Progressively higher-order terms in the expansion involve 
a larger number of flip-flops between the electron and 
nuclear bath spins. Consequently, higher-order terms are 
suppressed by the energy cost for such flip-flops, provided 
by the electron spin splitting b for an unpolarized nuclear 
bath. In particular, up to factors of order unity and an 
overall common prefactor, the size of the 2(n + l) th -order 
term is given by (see also Appendix A of Ref. 



£(2[n+l]) (t)oc 



I(I + 1)A 



(19) 



For a nuclear spin of order unity (/ ~ 1), the condition 
for the validity of a perturbative expansion in terms of 
Vg (i.e., the condition for convergence of the series in Eq. 
(fT8f ) is then given approximately by-£ ' 29 i 32 



b > A. 



(20) 



In the Born approximation, the self-energy E(t) is re- 
placed by the leading-order non- vanishing term in the ex- 
pansion of Eq. |[T8]): E(i) w t^(t). To understand the 
evolution of x t within the Born approximation, it is con- 
venient to introduce the function tp{t) = dt'T,^ 2 '(t'), 
which allows us to rewrite the equation of motion (Eq. 
CE5)) asi£ 



-i(V(0) 



Atj)x t + -R(t), 



(21) 



where R(t) = i J Q dt'ijj(t — t')xp. In a standard Born- 
Markov approximation, the dynamics induced by R(t) 
are neglected. The real part of ip(0) cancels any remain- 
ing precession: —Reip(Q) = Auj (see Eq. (fT4]) ) and the 
imaginary part of ijj(0) gives rise to a purely exponential 
decay of Xt with decay rate T = — Inr^(O). In a previous 
calculation, it has been shown that the decay rate for this 
system within a Born-Markov approximation vanishes: 
r = 0.' Within a Born approximation, all non-trivial 
dynamics in the rotating frame are therefore induced by 
the non-Markovian remainder term R(t). The remainder 
term R(t) has been investigated in detail previousl y 7 ' 25 ^ 6 
and leads to a partial decay of the coherence factor of 
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orderi \R(t)\ — O (4) 2 on a time scale ~ N/A, with 

long-time power-law tails. In the rest of this article, we 
include the effects of the Born approximation in induc- 
ing the Lamb shift (Eq. (H])), but neglect the < 0(1 /N) 
corrections due to R(t) in the perturbative regime. Ad- 
ditionally, we will include the fourth-order correction to 
the self-energy in the expansion of Eq. (fl"8|) . which we 
show induces a more dramatic decay than the Born ap- 
proximation, albeit at a longer time scale. 

We now approximate the self-energy by including all 
terms at second and fourth order in electron-nuclear spin 
flip-flops 



E(t) « S (2) (i) + £ (4) (i)- 



(22) 



Inserting Eq. ([2"2"| into Eq. (fT5|) we find, neglecting the 
dynamics with amplitude suppressed by ~ 1/N in the 
perturbative regime due to R(t): 




Xi 

Auj 



i \ dt't^\t-t')x v , (23) 



dt£ (2) (i). (24) 



-Re 



The integro-differential equation (Eq. ([23]) ) is diffi- 
cult to solve, in general. However, in terms of Laplace- 
transformed variables, this equation becomes an alge- 
braic equation, which can be solved directly. Introducing 
the Laplace transform of some function f(t), 



m- 

we rewrite Eqs. ([23]) and ([24]) as 

xq 



dte~ st f(t), Rc(s)>0, 



x(s) = 
Auj w 



« + iE( 4 )(«) 
-Re£ (2) (s = 0+). 



(25) 

(26) 
(27) 



We have calculated the self-energy E in Appendix [A] in- 
cluding terms up to fourth order in Vg. 

The dominant contributions to E(s) occur for \s\ <C 
\lo„\ in the rotating frame (high-frequency, s ~ iuj„ in 
the lab frame). We have expanded the second- and 
fourth-order self-energies in this limit, as described in 
Appendix A. Corrections to this expansion arc smaller 
than the retained contributions by a factor of A/Nui n <C 
1. For explicit calculation, it is useful to specialize to 
the case of a homonuclear system (where jk — 1 for 
all k) and a two-dimensional quantum dot with Gaus- 
sian envelope function, leading to coupling constants^^ 
Ak = (A/N)e~ k / N . Performing the continuum limit 
(J2k ~^ I ^0 ano - evaluating the relevant energy inte- 
grals for a uniformly polarized nuclear spin system leads 
to 



E< 4 )(s- iAuj) = ia[F+(s)J + (s) + F_(s)J_(s) - s\28) 
Auj ~ -E( 2 )(.s = 0+) = |(c + + c_)£4, (29) 



Figure 1: (Color online) Contour used to evaluate the 
Bromwich inversion integral (in the rotating frame defined by 
Eq. (|12[) ~1, The dynamics of the coherence factor xt are deter- 
mined by a single pole at so = — I/T2 and three branch cuts 
(see main text). The pole is offset from the branch cuts by 
the Lamb shift Auj and the excitation bandwidth (separation 
between branch points) is given by the size of the hyperfine 
coupling to a single nucleus at the center of the quantum dot, 
A/N. 



with 



24 



(30) 



and where c± = 1(1+1) — ((m(m±l)}} are the coefficients 
introduced in Ref. @ with ((•••)) indicating an average 
over the mixture of /^-eigenvalues, described by Eq. |7]). 
Additionally, we have introduced the functions 



F±(s) 
J ± (s) 



s ± i- 



A 
N 



log [ s ± i 



log(s). 



(31) 
(32) 



After inserting Eq. J25]) into Eq. (gBJ, we find that the 
Laplace-transformed coherence factor has three branch 
points, and if the principle branch is chosen for all branch 
cuts and at large electron Zeeman splitting (6 S> A), 
there is one pole at s = sq (see Fig. [1]). These non- 
analytic features determine the dynamics of the coher- 
ence factor (see below) . The equation-of-motion method 
adopted in Refs. I321I331 bears some similarity to the cur- 
rent approach. However, the excitation bandwidth (dis- 
tance between branch points) found in Refs. I321I331 is half 
that found here (A/2N rather than A/N), leading to a 
difference (by a factor of 2) for relevant decay time scales. 
We comment on other differences, below. 
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Since the rotating frame is chosen to give sq = — I/T2 
purely real, we find the general result 

x t = P e- t/T2 (35) 

i 

with Po given by Eq. (|4"5")) . below. The coherence factor is 
characterized by two terms: an exponential, which domi- 
nates in the perturbative regime (A < b, see Fig. [5]), and 
a sum of branch-cut integrals, which give rise to modu- 
lations of the decay envelope and a dominant long-time 
power-law decay (see Fig. [3} . 



Figure 2: (Color online) Comparison of the contribution from 
the pole at so in Fig. [T]with exponentially decaying residue 
(blue dashed line) with the full fourth-order result, obtained 
numerically (red solid line). We take the initial condition xq = 
1, assume an unpolarized nuclear spin system [with u n — b, 
c± = §I(J + 1), which follows from (<m 2 )) = /(/ + l)/3 
if all Zeeman levels have equal population], and have chosen 
I = 3/2 and A/b = 1/3. 
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Figure 3: (Color online) Long-time decay. At long times, 
the exponential decay envelope is modulated by branch-cut 
contributions. Parameters are as in Fig. [2] 



IV. SPIN DYNAMICS 

We find the time-dependent coherence factor by eval- 
uating the Bromwich inversion integral 



xt = lim / dse st x(s) 

7^0+ 2ni 



(33) 



which can be rewritten in terms of an integral over the 
closed contour Cb and branch-cut integrals /3j, j = 
0,+,- (see Fig. P: 



Res [e st x{s),s = s ] - ^ 0j(t). 



(34) 



In Eq. we have applied the residue theorem to write 
the integral over Cb in terms of a residue at the pole sq. 



A. Envelope modulations and long-time decay 

From direct asymptotic analysis of the branch-cut 
integrals, we find (3±(t) oc l/t 3 , while f3 (t) tx 
1/i 2 at long times. Since the pole contribution de- 
cays exponentially, the leading long-time asymptotics 
of Xt are thus given by /3o(t). Evaluating the pref- 
actor, we find the long-time limit (valid for t 3> 
max[l/Ao;, (6a/ Aw) In \NAuj/6aA\}): 



iAujt 



{2naA/N - i&ujf 



f 2 



(36) 



which gives the long-time behavior of the coherence fac- 
tor with initial condition xq = 1 (see Fig. 0]): 



where 



Re[x t ] 



C = 



C cos (Auit 



60 



(2waA/Ny + A 

( AloN\ 

—2 arctan 



\2naA ) 



(37) 

(38) 
(39) 



The modulations at a frequency Ato in Eq. (|37j) can 
be understood on physical grounds in the following way: 
The short-time dynamics of the electron spin are con- 
trolled by nuclear spins near the center of the dot, which 
are coupled most strongly. The effective precession fre- 
quency of the electron spin is therefore renormalized by 
the shift Ato due to virtual flip-flop processes with nuclei 
near the center for most of the decay envelope. The long- 
time decay, however, is controlled by weakly-coupled nu- 
clei far from the center of the dot, which cannot strongly 
shift the electron-spin precession frequency. The long- 
time dynamics therefore occur at the 'bare' precession 
frequency u> n . The difference in frequency between the 
dominant (short-time) and sub-dominant (long-time) be- 
havior leads to a relative beating at the frequency differ- 
ence Aoj. We note that the physical origin of this en- 
velope modulation is completely different from the more 
typical case of electron spin-echo envelope modulation 
(ESEEM), which is often observed in systems with an 
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Figure 4: (Color online) Comparison of the full numerical 
branch-cut integral — R e with the long-time asymp- 

totic expression given by Eq. (|37[) . Parameters are as in Fig. 

m 



anisotropic hyperfine interactioni 44 i 45 i 46 i 47 The modula- 
tions described here occur even in the present case of a 
purely isotropic interaction and without spin echoes. The 
two cases can be experimentally distinguished through a 
difference in the magnetic-field dependence of the mod- 
ulation frequency. Finally, we note that the modulations 
found here are reminiscent of modulations in branch-cut 
contributions that have been highlighted previously for 
the spin-boson model4£ 

Another striking feature of Eq. |37|) is the long-time 
power-law tail. This differs from the long-time expo- 
nential decay found by other author o 28 ! 30 ' 31 using re- 
summation and re-exponentiation techniques. The same 
long-time power-law decay (oc 1/t 2 ) has previously been 
predicted in Refs. I32ll33l based on an equation-of-motion 
method, but without mention of the phase shift or enve- 
lope modulations predicted by Eq. (l37|). To compare the 
results given here directly with Refs. 13 211331 . we consider 
the case of / = 1/2 with an unpolarizcd nuclear-spin 
bath. This gives 



C = 



(A/N) 2 
-tt + O 



O 
A 



(40) 
(41) 



Although the power law found here matches that re- 
ported in Refs. 32,33), and the modulations or phase 
shift can be ignored in the limit A/b <C 1, the prefactor 
C (which C = O [(b/A) 2 ] in Refs. [HlEl) is qual- 
itatively different. In particular, here we find that the 
power-law contribution with modulations can have sub- 
stantial weight (of order unity) in the perturbative regime 
A < b. This is clear from Fig. [U where we show that the 
branch-cut contributions can contribute approximately 
10% of the total decay amplitude. 



B. Decay shoulder 



For small t, we perform a Taylor-series expansion of 



x t - 



x t 



xo + x t + - x t + 



(42) 



From Eq. (|23[) and the initial value theorem we find 
Xo = and Xq = — iE^ 4 '(t = 0) = — tlinis^oo s& 4 \s). 
Inserting Eq. (j A18|) for £( 4 )(s) and choosing xq = 1 
gives: 



x t ^ 1 



t2> 



2ul 



(43) 



Eq. (l43l) gives the same short-time decay reported in 
Ref. [2f|, which was taken to describe a Gaussian coher- 
ence decay: x t ~ xq exp [— (t/T^) 2 ] • Here we note that 
the Gaussian approximation is only valid for times less 
than the actual decay time (t <§; r) in the perturbative 
regime (b > A) since the dominant decay is exponential 
in this regime for a typical (two-dimensional parabolic) 
quantum dot, as emphasized in Ref. 1291 . and illustrated 
here in Figs. [5] and 03 Re-exponentiation also fails for 
the fourth-order solution at lower magnetic field, as we 
show in Sec. [V] below. 

For a uniform unpolarized nuclear spin system, and 
for an electron with Gaussian envelope function in two 
dimensions, we find 



6^2 fb\ I N 



(44) 



We compare the initial decay found using this formula 
with the full non-Markovian solution in Fig. [5j While the 
short-time decay shoulder is well-described by a Gaus- 
sian, the full decay envelope is much better described by 
the dominant exponential solution (sec Fig. [5]). At larger 
Zeeman splitting b, the distinction between Gaussian and 
exponential becomes significantly more pronounced. 



C. Exponential decay 

We evaluate the residue at the pole Sq in Fig. Q] giving 

1 



Po 



(45) 



-1/T 2 



For a two-dimensional parabolic quantum dot, with 
an unpolarized nuclear system, we find Pq = 1 + 

O (J4] 2 In [4] ) . Thus, when A < 6, a Markov approx- 
imation is justified (resulting in a dominant exponen- 
tial decay), in agreement with the conclusions of Refs. 

In the Markovian regime, the decay rate for the ex- 
ponentially decaying pole can be determined through^ 
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Figure 5: (Color online) Decay shoulder. Here we show a 
comparison of the short-time quadratic decay with the full 
non-Markovian envelope for 1 — 3/2 and A/b — 1/3. 



I/T2 = — ImEW(0 + ). From the self-energy given in Eq. 
(|XT8|> . this gives 



1 

^2 



7TC + C_ 



(46) 



k.k' 



Here, the Markovian decay rate depends on the density of 
states for pair flips at an energy determined by the Lamb 
shift Aw oc A/u>„. The presence of Aw in the energy- 
conservation condition can be understood physically as 
arising from the rapid initialization that we assumed, giv- 
ing rise to the product-state initial condition (Eq. Q). 
At the instant the flip-flop interaction Vg is 'turned on', 
the electron spin experiences only the bare precession fre- 
quency w n . However, after some interaction time scale, 
the renormalized precession frequency w n + Auj gives 
the electron energy splitting, and so the correct energy- 
conservation condition for nuclear-spin pair flips contains 
the difference of these two quantities (i.e., Alu). The de- 
pendence on Aw shown in Eq. (|46[) results in an interest- 
ing (in general, non-monotonic) dependence of I/T2 on 
magnetic fields In contrast, the effective-Hamiltonian 
approach that was adopted in Rcfs. 27.28.29.30.31 incor- 
porates the Lamb shift as an additive constant directly 
into the electron-spin Zeeman splitting, and so it does 
not enter into the formula for I/T2. The Lamb shift that 
comes out of the same procedure used here, but start- 
ing from the effective Hamiltonian^ has a dependence 
Aw oc (A/u) n ) 2 1 so the effective-Hamiltonian treatment 
does not recover the correct magnetic-field dependence 
given by Eq. (|4"6")) . although the leading-order Aw ~ 
behavior is recovered correctly. This is not surprising, 
since the effective Hamiltonian is only strictly valid to 
leading order in A/u> n . Through explicit calculation at 
higher orders, we have checked that the leading correc- 
tion to the Markovian decay rate at sixth order in Vg 

is O (A/ojn) 4 (see Appendix IB]) , and so the expression 

given here is at least correct up to and including terms 

of order O \{A/ui n f 



Figure 6: (Color online) Decoherence rates F; from Eq. 
(|49p and total decoherence rate I/T2 = Li for an elec- 
tron spin in a GaAs quantum dot containing N = 10 nu- 
clei with <?-factor \g\ = 0.4. The decoherence rate shows a 
non-monotonic behavior, reaching a pronounced maximum. 
This is in contrast to the leading-order result,— and in con- 
trast to the results from other higher-order expansions ~ 
We have used hyperhne coupling constants A Ga = 74 /xeV, 
A Ga — 96^eV, A " Aa = 86 fieV and relative abundances 
i/69 Ga = 0.3, V71 g„ = 0.2, V75a 3 = 0.5, which have been esti- 
mated in Ref. Eg (see Table 1 of Ref. (H). 



Specializing to a two-dimensional quantum dot with 
Gaussian envelope function and evaluating the energy 
integrals in the continuum limit gives 

T = ^V'v W + 2e2 ) ( 1 - £ )' ( 47 ) 

T 2 3(c + + c^y N 

where 

c+ + c_ 



A_ 



(48) 



First, we note that the sub- leading contribution to I/T2 
in Eq. (|47| (oc e 3 ) is suppressed only by one power of 
A/ui n (up to corrections of order unity). Second, this 
sub-leading correction has the opposite sign of the lead- 
ing (~ e 2 ) term, potentially leading to a non-monotonic 
dependence of I/T2 on the electron Zeeman splitting 
when e ~ 1. This non-monotonic dependence can be 
understood in the following way: As the electron Zee- 
man energy decreases from a large value b 3> A, the 
perturbative Lamb shift Alu oc 1/6 increases, eventually 
reaching the edge of the band of single nuclear pair-flip 
excitations when Aw ~ A/N , at which point there are 
no more energy-conserving flip-flop processes. For still 
lower magnetic fields, higher-order processes are required 
to conserve energy, but we find that these processes are 
further suppressed by the small parameter c+c_ for a po- 
larized nuclear spin system (c+c_ oc (1 — p 2 ) for nuclear 
spin / = 1/2) (see Appendix [B]) . The qualitative non- 
monotonic magnetic-field dependence described by Eq. 
(|4"7|) will therefore apply at least in the case of a polar- 
ized nuclear spin system, even when higher-order terms 
in Vg are taken into account. 
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Figure 7: (Color online) Total decoherence rate I/T2 = X]j Ti 
(solid line) with Ti from Eq. (|49|l for an electron spin in a 
In^Gai-j; As quantum dot containing N = 10 nuclei with g- 
factor \g\ = 0.5 and In doping x = 0.3. We show individual 
contributions from F75 As (dotted) and Ti n (dashed). Contri- 
butions from the gallium isotopes are not visible on this scale. 
For this plot, we have taken hyperfine coupling constants for 
the gallium and arsenic isotopes as in Fig. [6] and have used 



A = A 
of Ref. [H). 



A In = HO^eV from Ref. \M (see Table 1 



It is straightforward to extend the analysis of this sec- 
tion to the case of a heteronuclear system. Provided the 
difference in nuclear Zeeman energies exceeds the exci- 
tation bandwidth (|(7j — jj)b\ > A/N), the decoherence 
rate is given by a sum of contributions from each nuclear 
species i: I/T2 = X^^i, where, for a two-dimensional 
quantum dot with Gaussian envelope function, we find 



irisfai (ef - 3e, 



A' 1 

2) —9(1 



NAlu 



Alo 



E< 



- £i 



A % A 1 



8uj„ 



(49) 
(50) 



Figure 8: (Color online) Decoherence rate I/T2 for an electron 
spin in a In^Gai-a; As quantum dot containing N = 10 nuclei 
with (/-factor \g\ = 0.5 and In doping of x — 0, x — 0.3, and 
x — 0.5. Hyperfine coupling constants are as given in the 
caption of Fig. [7] 



V. NON-PERTURBATIVE REGIME: b < A 

Although the expression we have given for the self- 
energy is strictly valid only in the perturbative regime 
(b 3> A), here we explore the non-Markovian dynamics 
of this solution outside of the regime of strict validity and 
comment on where the results become unphysical. 

We find the positions of poles and evaluate residues 
and branch-cut integrals numerically to find the coher- 
ence factor in this regime. We consider the case of an 
unpolarizcd homonuclear spin system with spin / = 3/2, 
appropriate to GaAs. As the electron Zeeman splitting 
b is lowered from b ^> A, we find there is a critical value 
of b (near b ~ 2A), below which there is a second pole 
(at s = si) with exponentially-decaying residue. The co- 
herence factor xt is then given by a sum over two pole 
contributions and three branch-cut integrals: 



and where we have introduced 



= E E &(*)• 



(52) 



.7=0,+,- 



"■-^(s)'- (51) 

with coefficients c l ± = P(P + 1) - ((m 4 (m l ± 1))) for 
each isotopic species i. We show the magnetic- field de- 
pendence of the decoherence rate 1/T 2 from Eq. (|4"9"|) in 
Fig. [6] for a GaAs quantum dot, and in Fig. [7| for an 
InGaAs quantum dot with a typical indium doping of 
x = 0.3. The dependence of the I/T2 curve on indium 
doping x for an InGaAs quantum dot is illustrated in 
Fig. [51 where we see that the position of the maximum 
in the I/T2 curve depends strongly on the concentration 
of the large-spin isotope (indium) . An experimental con- 
firmation of this dependence of the maximum in 1 /T 2 as 
a function of indium doping would be a strong confirma- 
tion of this theory. 



For b = 2 A (top panel of Fig. 0, there are two 
exponentially-decaying pole contributions, giving rise to 
a bi-exponcnatial decay with strong envelope modula- 
tions corresponding to the difference in the imaginary 
part of the two poles. At smaller Zeeman energy A/2 < 
b < A (e.g. b = A in the center panel of Fig. [5]), the 
pole at s = so leaves the continuum band and merges 
with the imaginary axis, leading to a constant contribu- 
tion Po(t) = Pq, independent of t. For still lower Zeeman 
splitting b < A/2 (b = A/2 in the lower panel of Fig. O, 
the second pole at s = s% leaves the continuum band at 
lower frequency and also merges with the imaginary axis. 
In this regime, the only decay in the fourth-order solu- 
tion is due to the small contribution from branch cuts, 
although envelope modulations remain. 

The effects in the two lower panels of Fig. [S] demon- 
strate the danger of re-exponentiation of short-time be- 



9 




Re [Ei -Pi(t) - Ej Bj(t)] 
Re [^ = - t/T 2] 



10 

i (iV/A) 



20 



Figure 9: (Color online) Decay envelopes calculated from nu- 
merical evaluation of branch-cut integrals and pole contri- 
butions for a two-dimensional quantum dot with Gaussian 
envelope function (solid lines). Dashed lines show the contri- 
butions from the dominant pole at s = so . When the electron 
Zeeman splitting b is below some critical value b < b c ~ A, 
a second exponentially-decaying pole appears, leading to a 
biexponential decay with strong envelope modulations (top, 
b = 2A). When b is decreased further, the dominant pole 
moves to the real axis (middle, b = A). At still smaller values 
of b, the sub-dominant exponential pole has a vanishing de- 
cay rate (bottom, b = A/2), leading to sustained oscillations. 
All plots correspond to an unpolarized narrowed nuclear bath 
with / = 3/2. Insets illustrate the approximate relative po- 
sitions of poles (circles) and branch points (crosses) in each 
case. 



havior for a system where strong non-Markovian (history- 
dependent) effects become important. The non-decaying 
fractions shown in Fig. [5] are, however, unphysical conse- 
quences of extending the solution to a regime of electron 
Zeeman splitting where it does not apply. We expect 
higher-order corrections to the self-energy to broaden the 
continuum band as higher-order nuclear pair flips are in- 
cluded, resulting in several exponential decay time scales 
as the electron Zeeman energy is lowered. Nevertheless, 



we have found that processes that can broaden the con- 
tinuum band will be suppressed even at small electron 
Zeeman splitting b < A, provided the nuclear-spin sys- 
tem is polarized (see Appendix [Bj , and so some of this 
behavior will survive at least for a polarized nuclear-spin 
environment. Whether perturbation theory can be con- 
trolled at any polarization for b < A through an ad- 
equate resummation of relevant terms, as suggested in 
Refs. l3Cill3ll . is still unclear with the present method. 



VI. CONCLUSIONS 

We have investigated transverse-spin dynamics for an 
electron confined to a quantum dot, interacting with a 
bath of nuclear spins via the Fermi contact hyperfine 
interaction. Using one unified technique, we have re- 
covered results that have previously been reported using 
several different methods. These results include an initial 
partial decay, followed by a quadratic shoulder, a domi- 
nant exponential decay, and a long-time power-law tail. 
Our results for the long-time behavior differ from those 
of Refs. H^HlllSIIl Here, we have found a long-time 
power-law decay (~ 1/i 2 ), in contrast to the long-time 
exponential decay found by those authors. While the de- 
cay law ~ l/t 2 matches that found previously using an 
equation-of-motion approach ) 32 i 33 the prefactor found in 
the present work has a qualitatively different dependence 
on magnetic field. In contrast to earlier works, which ar- 
gue in favor of a regime of Gaussian decay ; 27 i 28 here we 
find that re-exponentiation of the short-time quadratic 
decay shoulder is never justified. In the perturbative 
regime b > A, the system is Markovian^ being well- 
described by a single-exponential decay. As the electron 
Zeeman splitting is lowered to b < A, we find strong non- 
Markovian effects (sustained oscillations and multiple de- 
cay rates), which once again invalidate re-exponentiation 
of the short-time decay shoulder. 

In addition to recovering previous results, we have 
found qualitatively new behavior, including modulations 
of the decay envelope and sub-leading corrections to 
the decoherence rate for the dominant exponential de- 
cay. Our calculation gives an interesting non-monotic 
dependence of the decoherence rate I/T2 on magnetic 
field. These two results (envelope modulations and a 
non-monotonic dependence of the decoherence rate on 
magnetic field, both of which should be readily accessi- 
ble in experiment) are not recovered in dynamics under 
a leading-order effective Hamiltonian, suggesting caution 
should be exercised in interpreting results of high-order 
expansions involving the effective Hamiltonian. 
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Appendix A: SELF-ENERGY EXPANSION 

Here we give the explicit self-energy up to fourth or- 
der in Vff. The full self-energy supcropcrator is given by 
E(s) = ZW(s) + Z( 4 >(s) + 0(V£), where 



G(s) 



-iTr(S + L v G(s)L v S- Pl (0)), 
1 



s + iL a ' 
l O = [H ,O], 
L v O = [Vs,0], 



(Al) 

(A2) 

(A3) 
(A4) 



and the fourth-order result is 

S( 4 )( s ) = iTr [S+ (1 - iL QG(s)) L v G(s)L v Q 

xG(s)L v G(s)L v S-pi(0)]. (A5) 

More explicitly, using Eq. ([5]) for the initial nuclear state 
we find SW(s) = J2i Pu^f\s); where 



zf( s) 



4e 



,+ 



km 

l h+ } ni k l h ~} kni 
s + i5uj kn + iuL 



(A6) 



kin 



Ef( S ) 



16 E [ /l -]feife 2 [ /l +]fe 2 fc 3 [ ;i -]fe 3 «i [ G T]fe ini [ G +}k 2 n l i G ^k 3ni i 1 - S nik 2 ) 



kik 2 k 3 



+ [h-]n,ki \- h +]k 1 k 2 i h ~]k 2 k 3 i h +}k 3ni i G l} ni k 3 \- G +\ ni k 2 [ G l] nikl C 1 ~ <W 2 ) 

+ [fc-U [*+U, [M fc2fc3 M*,^ (l - « (^o + L 2fc2 ~ [^o + ] ni „ ; ) [<?+]*.*.) 
J^U [G+U + [G T ] fc3 , 2 [G + ]„ ;fc2 [GJ^ J (1 - 5 nife2 ) 

+ (j G T] fc3fc2 + [ G lk 2fel ) [ G -] fe3fel (t^U + [GlUj] } ■ ( A? ) 



Here, we denote matrix elements [/i ± ] r 



Further, (5w„ fc = \ (/i* - hjj!) and w^ fc = (n| |n) - 
a; 7 |fc), where to 1 = b^2 k ^kl^i an d we have introduced 



[ G a fcfc' — 



(A8) 



with 



ioliXTIIA) 

ioit> aiifc) (fc'i 

io|T) (Tllfc) (A'| 

Lo\i)(i\\k) (fc'| 



[L+] kk ,\i)(^\\k)(k'\, (A9) 
[io] fefc Jt>ai|fc)(fc'|, (A10) 
|T)(T||fc)(fc'|, (All) 



L 



T I 
^0 



kk' 



kk' 



u)aiifc)(fc'|. (Ai2) 



The dominant contributions to the self-energy occur 
at high frequency (s ~ i^> n ) m the lab frame. For 

\s — iuj n \ <C u> n , and uj n 3> 5uj nk , uj„ 3> ui^ k , we have: 



\G 



[G 



l\kr. 



1 



1 + 



.4 



Nuj n 



(A13) 



which allows us to approximate Eqs. (|A6[) . (| A7[) by their 
high-frequency forms for a uniformly polarized system. 
We additionally go to the rotating frame; from the defi- 
nition of £ in Eq. lfl~3"l) , we have 



S(s + iuj n ) = £(s — iAw), 



(A14) 



which gives 



rE 1 



4w„ 



(A15) 

In the above expression, the sum over i indicates a sum 
over different nuclear-spin isotopes with abundances Vi 
and hyperfine coupling constants A\ . The high-frequency 
form of the fourth-order self-energy is then: 
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£( 4 )(s-iAcj) = 




-12 



IX-,' 



»7y 



1 



(A16) 



where x^ 2 = (A J fci - AjJJ/2, 7 y = 6(7, - 7,-), and the 



coefficients Cj_ are: 

<4=J 4 (Ji + l)-«m(m±l))) 



(A17) 



With the average ((■■■)) defined in Eq. 0. For a 

ij 

-12 — x 12 — 

1. In this case, 



homonuclear system, we have 7^ = and x 



(A kl - A k2 )/2, and replace 

assuming a uniformly-polarized nuclear-spin system, the 
self-energy is given simply by 



S (4) (s-iAw) 



■ C+C- 

' 4wi ^ s 

n k,k' 



i{A k 



Ay. 



(A18) 



This self-energy differs from that found previously at 
leading order in an effective- Hamiltonian treatment,— 
where the Lamb shift Atu is incorporated directly into the 
bare precession frequency w n . In addition, we stress that 
the more general self-energy for a heteronuclear system 
(Eq. (|A16|) ) is not recovered with the effective Hamilto- 
nian (compare with Eq. (C19) of Ref. l29h . 

Assuming a two-dimensional parabolic quantum dot 
(with Gaussian envelope function) leads to coupling con- 
stants A k = (A/N)e k/N (see, e.g., Ref. [!§). Performing 
the continuum limit, i.e., replacing ^2 k k — ► J dkidk^, 
and evaluating the resulting energy integrals, we arrive 
at Eq. (|28|) of the main text. 



Appendix B: HIGHER-ORDER CORRECTIONS 

All results in this article are valid up to fourth order 
in the electron- nuclear flip-flop terms Vg. As the elec- 
tron Zeeman splitting is lowered from b A, higher- 
order corrections to the self-energy may become rele- 
vant. In this Appendix, we give explicit conditions under 
which higher-order corrections may be neglected, even 
for b ~ A. As in Appendix \^ the self-energy at any 
order may be approximated by its high-frequency form 
(at s ~ iu) n ) whenever A/Nb <C 1. This allows for a 
significant simplificaiton in the high-order expansion in 
terms of Vg. With corrections to the self-energy that are 
smaller by factors of order 1/N <C 1 and A/Nb <C 1, 
we find the high-frequency form of the self-energy to be 
given by 



E(s) ~ — iTii 



G7 Q 



1a-1 



2 u 



1 



-P/(0) 



, (Bl) 



where we have introduced 

1 



G+ = 



s - -L H 

6 2 Lu 



o = -t^-G+ (Hi + H fl ) 

4w n 



(B2) 
(B3) 



defined in terms of the superoperators (which act on an 
arbitrary operator O): 



H lO = h+h-O, 
H r O = Oh-h+, 
L+O = {u>,0}, 



(B4) 
(B5) 
(B6) 



where {,} indicates an anticommutator. 

The high-frequency form of the self-energy can now 
be found directly from Eq. (|Bip with a more moderate 
constraint on the electron Zeeman splitting (A/Nb <C 1). 
A direct evaluation of Eq. (|B1[) at arbitrary order and 
resummation is non-trivial, but we can generate arbitrary 
higher-order terms with the geometric series: 



1 



1 



= l-<7 



(B7) 



Every factor of a is associated with a nuclear-spin pair 
flip, giving rise to a factor of c + or c_, which depend on 
the nuclear polarization, and a factor of the small pa- 
rameter A/ui n . The term at (2fc) th order in Vg contains 
k factors of a, and consequently k powers of A/ui n . This 
suggests that the sixth-order self-energy can in general 
give corrections of order ~ (A/uj n ) 3 , which may mod- 
ify the sub-leading corrections of this size given by the 
Markovian decay formula (Eq. (fi5|) ). However, by direct 
calculation using the above expansion, we find the lead- 
ing contributions to the Markovian decay rate at sixth 
order to be 



ImS (6) (s 



O 



(B8) 



Furthermore, we find that the E^ 6 ) corrections do not 
lead to a broadening of the continuum band. The first 
non-vanishing corrections to the Markovian decay rate 
which do lead to a broadening of the continuum band 
contain two nuclcar-spin pair-flip excitations. These 
terms occur first at eighth order in Vg and are sup- 
pressed by the factor (c+c_ ) 2 -which is smaller than the 
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fourth-order corrections by a factor c+c_ for a polarized 
nuclear-spin system (e.g. c + c_ oc (1— p 2 ) for nuclear spin 
/ = 1/2). This result demonstrates that the qualitative 
decrease in the decoherence rate at low electron Zeeman 
splitting shown in Figs. [6l [3 and [8] will not be signif- 



icantly modified by higher-order corrections, at least in 
the case of a large polarization, where perturbation the- 
ory still applies at a smaller value of the electron Zeeman 
splitting. 
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